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Abstract 

The latest results from atmospheric and accelerator neutrino experiments indicate that the normal 
neutrino mass ordering rrii < m 2 < a maximal leptonic CP-violating phase S = 270° and the second 
octant of neutrino mixing angle ^23 > 45° are favored. In light of new experimental results, we update 
previous phenomenological studies on two-zero textures of the Majorana neutrino mass matrix Af^, in 
the flavor basis where the charged-lepton mass matrix is diagonal. When the Icr ranges of neutrino 
mixing parameters are taken into account, only four (i.e., 2 and B2 4) among seven two-zero patterns 

of show the aforementioned features of neutrino mass spectrum, mixing angle ^23 and CP-violating 
phase (5, and thus are compatible with the latest neutrino oscillation data. The correlative relations 
among neutrino masses and mixing parameters have been derived analytically for these four patterns, 
and the allowed regions of neutrino mixing angles and the CP-violating phase are also given. Possible 
realizations of four viable two-zero textures via non-Abelian discrete flavor symmetries are discussed. 
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1 Introduction 


In 2011, the accelerator neutrino experiments T2K [T] and MINOS [2] observed neutrino oscillations in 
the —)■ t'g appearance channel, indicating a nonzero 9^^ at the 2a significance level, which is also 

consistent with the result from the reactor neutrino experiment Double Chooz [3]. In 2012, the Daya Bay 
collaboration [3] observed the disappearance of from nuclear reactors and discovered a relatively large 
value 0^3 ~ 9° with a high statistical signihcance, which was further confirmed by the RENO [5] and 
Double Chooz [ 6 ] experiments. Thanks to three leading reactor neutrino experiments, Daya Bay, Double 
Chooz and RENO, the value of the smallest mixing angle 6i^ has been measured more precisely than 
those of the other two large mixing angles 023 - The remaining important issues in experimental 

neutrino physics are to determine whether neutrino mass ordering is normal (i.e., < m 2 < m^, denoted 

as NO) or inverted (i.e., m 3 < < m 2 , denoted as 10), to discover leptonic CP violation and measure 

the CP-violating phase 5, and to pin down the octant of @23 (i-®-) ^23 < ^23 > 45°). Any important 

experimental progress will shed light on the underlying flavor structure of lepton mass matrices, and thus 
help us explore the origin of neutrino masses and flavor mixing. 

Very recently, the T2K collaboration has published a combined analysis of the data in both —>■ 

disappearance and —>■ z/g appearance channels [7]. If the reactor neutrino data relevant for 6^^ 
are included, the likelihood maximum is reached at 5 = 270° and sin^ ^23 = 0.528 (i.e., @23 = 46.6°) 
in the NO case. In addition, the latest preliminary results of appearance from the NOzzA 

experiment also favor NO at the 2cj level, and point to <5 ~ 270° if @23 = 45° is fixed [ 8 ]. Apart from 
the long-baseline accelerator neutrino experiments T2K and NOzzA, the atmospheric neutrino experiment 
Super-Kamiokande has presented the latest results on neutrino oscillation parameters [9], showing a weak 
preference for NO, a nearly maximal CP-violating phase b = 240°, and the second octant @23 > 45°. 
The preference for NO arises from an excess of up-going zZg-like events, which could be interpreted by 
the resonant flavor conversion associated with ^33 in the neutrino channel. The combined analysis of 
accelerator and reactor neutrino data has been performed in Ref. m, where the values of 5 around 90° 
are found to be disfavored at between 2cr and 3a, depending on neutrino mass ordering and the octant of 
622 - Similar conclusions are also reached in Ref. [11]. Although the statistical significance of the above 
results is below 3a for each experiment, we have the same indication from both accelerator and atmospheric 
neutrino oscillation experiments, for which the neutrino sources, oscillation physics and detection methods 
are completely different. Therefore, it is very interesting to study the implications of a combination of 
NO, (5 = 270° and @23 > 45° for the flavor structure of the Majorana neutrino mass matrix^ Based on 
the possible flavor structure, one can go further to explore the dynamical origin of neutrino masses and 
underlying flavor symmetries. 

In this paper, we update the phenomenological studies of two-zero textures of Majorana neutrino mass 
matrix in Ref. [12]. The motivation of such an update is two-fold. Eirst, in the flavor basis where the 
flavor eigenstates of charged leptons coincide with their mass eigenstates, it is straightforward to show that 
two-zero textures of Majorana neutrino mass matrix contain only five real parameters, which should be 
confronted with nine observables, including three neutrino masses, three flavor mixing angles, and three 
CP-violating phases. As a consequence, they are very predictive, and can readily be tested in neutrino 
oscillation experiments. Second, Ref. [ 12 ] and others [13 [H] were inspired by the discovery of a nonzero 
0 ^ 3 , so they did not focus on the mass ordering, a maximal CP-violating phase or the octant of 023- 

^In the present work, neutrinos are assumed to be Majorana particles, as they indeed are in a class of seesaw models for 
neutrino mass generation. 
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particular, only the parameter space around 5 = 90° or satisfiying ^23 < 45° is concentrated on in Ref. [12], 
which is now inappropriate in light of recent T2K, Super-Kamiokande and NOi/A results. Furthermore, 
we have derived a novel analytical formula for the correlation among neutrino mixing parameters, which is 
applicable to a special set of two-zero textures and turns out to be much simpler than that existing in the 
literature. Among seven two-zero textures, we find that only four (i.e., 2 B 2 4 ) remain consistent 

with the latest neutrino oscillation data if the la ranges of mixing parameters are taken into account. 

The rest of this paper is structured as follows. In Sec. 2, the important and basic formulas for two- 
zero textures are summarized. Sec. 3 is devoted to analytical and approximate results of theoretical 
predictions, and a complete numerical analysis of the allowed parameter space. We have found a novel 
and useful formula that has not yet been noticed in previous works. It turns out that only four out of 
seven patterns (i.e., 2 and B 24 ) are compatible with the latest neutrino oscillation data. The precise 

determination of S and 6*23 could further single out a unique pattern of two-zero textures, or exclude all of 
them eventually. In Sec. 4, we present some examples to illustrate how to realize those two-zero textures 
through non-Abelian discrete flavor symmetries. Finally, our main conclusions are summarized in Sec. 5. 

2 Two-zero Textures 

Since neutrino masses cannot be accommodated within the standard model of elementary particles (SM), 
one can regard the SM as an effective theory at the electroweak scale and introduce higher-dimensional 
operators, which fully comply with the SM gauge symmetries. As already pointed out by Weinberg m. 
the lowest-dimensional operator composed of lepton and Higgs doublets is of dimension five, 

which is unique and leads to lepton number violation by two units. After spontaneous gauge symmetry 
breaking, neutrinos acquire a tiny Majorana mass term from this five-dimensional operator. Such an 
operator can actually be realized in renormalizable models, e.g., the seesaw models, in which three singlet 
right-handed neutrinos [El dZl [El [El EQ] or one triplet scalar [211 \22l EHl EH EH ES] or three triplet 
fermions [27| are introduced and the mass scale of these new particles is extremely high. 

In the seesaw framework, the mass terms relevant for lepton mass spectra and flavor mixing at low- 
energy scales can be written as 

~ 2 ’^L^u’^L + h-C. , (I) 

where and stand for the charged-lepton and Majorana neutrino mass matrix, respectively. In the 
flavor basis, where = diag {m^, m^, is diagonal, the lepton flavor mixing matrix V stems from the 
diagonalization of a complex and symmetric matrix M^, namely, V'^M^V* = = diag{m^,m 2 ,nig}. 

Here m„ (for a = e,fi,T) and (for i = 1,2,3) denote respectively the charged-lepton and neutrino 
masses. For later convenience, we decompose the mixing matrix as H = U ■ P, where P = diag{e'^, 6'°^, 1} 
with {p, a} being two Majorana CP-violating phases, and U can be parametrized in terms of three mixing 
angles {^ 125 ^ 135 ^ 23 } Dirac CP-violating phase 5, namely. 



' C12C43 

^ 12^13 

5436 


u = 

“^12^23 “ ^125135236''^ 

+^12^23 “ ®12'5i3'®23®'^ 

^ 13^23 

(2) 


V+Si2'523 “ ^125130236^'^ 

~‘^ 12^23 ~ ^ 12 ^ 13 ^ 23 ^^^ 

^13 *^23 ) 



with Sjj = sin and = cos for ij = 12,13,23. The three mixing angles { 012 , ^ 13 , ^ 23)5 neutrino 
mass-squared differences Am^i = m^ — rnl and Amgi = m 3 — mf, and the Dirac CP-violating phase 6 can 
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be determined from neutrino oscillation experiments. In addition, the constraints on the absolute neutrino 
masses can be obtained from both beta-decay and neutrinoless double-beta decay experiments, while the 
latter could also provide us with some useful information on the two Majorana CP-violating phases. The 
latest observations of cosmic microwave background by the Planck collaboration have placed a restrictive 
bound on the sum of neutrino masses S = nii + m 2 + m^ < 0.23 eV at the 95% confidence level |28j . 

In this context, any useful experimental information on the mixing matrix V and neutrino masses 
(for z = 1, 2, 3) can be implemented to reconstruct the Majorana neutrino mass matrix via 



/ m3 
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0 ^ 
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where A;^ = A 2 = 771 , 26 ^^'^, and A 3 = m 3 . On the other hand, any constraints on M^, arising from 

flavor symmetries or empirical assumptions will result in specific predictions for neutrino mass spectrum 
and mixing parameters, which are ready to be tested in neutrino oscillation experiments. Now we follow 
the second way and assume two texture zeros in M^, which were first considered by Frampton, Glashow 
and Marfatia in Ref. [29] and thoroughly studied by Xing in a series of papers [301 ED [32], and by some 
other authors |33[ 1341 (35] . The two-zero textures have been reexamined after the experimental discovery 
of a nonzero ^33 in Refs. [laiiaiiiiEgEaEsiES]. 

The predictions from two-zero textures and their experimental tests can be found in previous works, 
however, we summarize the important formulas in this section and establish our notations. Since is 
symmetric, there are totally six independent matrix elements. If two of them are assumed to be zero, i.e., 
{M^)ab = = 0 for ab / cd, we obtain 


^U,lUbl + ^U,2Ub2 = -UasU,,, 

^U,iUdi + ^U,2Ud2 = -UM,. (4) 

^3 ^3 

From Eq. (4), one can immediately extract the neutrino mass ratios |30l ETl (321 [12] 


m 3 

m 2 ^ 

m 3 

and two Majorana CP-violating phases 

1 

P = 2 

1 

cr = - arg 


C = 


Ua,U,,U,2U,2-U,2Ub2Uc,Ud, 


Ua2Ub2UclUd^ - U,lU,,U^2Ud2 


Ua2Ub2UclUdl - U,,U,,U^2Ud2 


Ua,UbzUc2Ud2 - Ua2Ub2UcUd3 
Ua2Ub2UclUdl-UalU,^U,2U,2\ ’ 
UglUblUczUd, - 
Ua2Ub2UclUdl - U,lUblU,2Ud2^ 


( 5 ) 


( 6 ) 


Now it is instructive to count the number of free parameters in this scenario. The physical parameters 
include three neutrino masses {m 3 ,m 2 ,m 3 }, three mixing angles {^ 121 ^ 135 ^ 23 } three CP-violating 
phases {5, p, a}. Given the four constraints or relations in Eqs. (5) and ( 6 ), we are left with five degrees of 
freedom, which can be fixed by two neutrino mass-squared differences Am|^ and Am^^, and three mixing 
angles {^121 ^131 ^ 23 }- Therefore, the absolute neutrino mass (e.g., m 3 ) and CP-violating phases (d,/?, cj} 
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can be determined by the available neutrino oscillation data on Am|^} and {^ 12 )^ 13 )^ 23 }- To be 

explicit, neutrino masses can be found by using 


m 3 = 


Am 2 ^ 

^ 2_^2 


^2 = C”^3 , rrii = ^mg 


(7) 


and the Dirac CP-violating phase 5 can be fixed by 

_ Am2g 

"-|Amig| |l-e2| 


( 8 ) 


where the dependence of C and ^ on (i is implied in Eq. (5). The approximate analytical results of {5, p, cr} 
can be found in Refs. [301 [311 ESI Eg. In Ref. [ig, the global-fit results of neutrino oscillation data at 
that time have been used to check the viability of two-zero textures of It turns out that the following 
seven patterns survive current neutrino oscillation data: 





) A 2 





, B 2 


/a 0 

( c d 
\b d Oy 


(9) 


( 10 ) 


where the nonzero elements are in general complex. However, it is always possible to make three nonzero 
matrix elements real by redefining the phases of three charged-lepton fields. In this case, we have five free 
model parameters (e.g., three real elements {a, b, c} and one complex element d), as it should be consistent 
with our parameter counting. 

Motivated by the recent results from T2K, Super-Kamiokande and NOz^A, we reexamine those seven 
patterns by focusing on whether NO, 6 = 270° and ^23 > 45° are allowed simultaneously. If this is the case 
for a given pattern, we further completely explore the viable regions of neutrino masses and flavor mixing 
parameters. Unfortunately, the latest global-fit analysis of neutrino oscillation data has not yet included 
the preliminary results from T2K, Super-Kamiokande and NOz^A, although partial joint analyses of reactor 
and long-baseline accelerator neutrino data have been performed in Refs. mm For comparison, we 
collect the global-ht data from Ref. [4U] in Table 1. Due to the lack of a full statistical analysis, we quote 
the la intervals of Am^g and 0^2 from Table I in the NO case, while the 68 % Bayesian credible regions 
of Am 32 , ^ 23 ) ^13 S are taken from Ref. [g, in which a Bayesian analysis of both disappearance and 
appearance data from T2K and reactor neutrino data has been performed. The oscillation data used in 
our numerical calculations are summarized as follows 


32.73° < 6»i2 < 34.26° , 8.61° < 0^3 < 9.56° , 44.43° < 023 < 49.78° , 202° < <5 < 338° (11) 

and 

7.33 X 10“^ eV^ < < 7.69 x 10"^ eV^ , 

2.47 X 10“^ eV^ < Amii < 2.69 x 10"^ eV^ . (12) 

The best-fit values of three mixing angles (0i2 = 33°, 0^3 = 9°, 023 = 47°), two mass-squared differences 
(Am^g = 7.50 X 10“^ eV^, Am^^ = 2.58 x 10“^ eV^), and the CP-violating phase 6 = 270° will also be 
used for illustration. 
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Table 1: The best-fit values, together with the Icr, 2a and 3cr intervals, for three nentrino mixing angles 
{^ 12 )^ 13 )^ 23 }) two mass-squared differences or Am| 2 } and the Dirac CP-violating phase 6 

from a global analysis of current experimental data m- Two independent global-fit analyses can be found 
in Refs. HD 02], which are in perfect agreement with the resnlts presented here at the 3a level. 


Parameter 

Best ht 

Icr range 

2a range 

3fj range 

Normal neutrino mass ordering (m^ < m 2 < m^) 

012 /° 

33.48 

32.73 — 34.26 

31.98 — 35.04 

31.29 — 35.91 

013/° 

8.50 

8.29 — 8.70 

8.08 — 8.90 

7.85 — 9.10 

023/° 

42.3 

40.7 — 45.3 

39.1 — 48.3 

38.2 — 53.3 

v° 

306 

236 — 345 

0 ~ 24 © 166 — 360 

0 — 360 

Amii/[10-® eV^] 

7.50 

7.33 — 7.69 

7.16 — 7.88 

7.02 — 8.09 

Am|;^/[10“^ eV^] 

+2A57 

+2A10 — +2.50A 

+2.363 — +2.551 

+2.317 ~ +2.607 

Inverter 

. nentrino mass ordering (rug < m^ < m 2 ) 

012 /° 

33.48 

32.73 — 34.26 

31.98 — 35.04 

31.29 — 35.91 

CO 

0 

8.51 

8.30 — 8.71 

8.09 — 8.91 

7.87 — 9.11 

023/° 

49.5 

47.3 — 51.0 

45.1 — 52.5 

38.6 ~ 53.3 

v° 

254 

192 — 317 

0 — 20 © 130 — 360 

0 — 360 

Amii/[10-® eV^] 

7.50 

7.33 — 7.69 

7.16 — 7.88 

7.02 — 8.09 

Am|2/[10“^ eV^] 

-2.449 

-2.496 -2.401 

-2.543 -2.355 

-2.590 -2.307 


3 Analytical and Numerical Results 


Before proceeding with a numerical analysis, we briefly comment on the analytical results presented in 
Refs. [301 [311 [ 32 I [ 12 ] , where the standard parametrization of U in Eq. (2) has been inserted into Eqs. (5) 
and (6) to derive the neutrino mass ratios and CP-violating phases in terms of mixing angles. Of course, 
this is a correct and straightforward way to achieve an analytical understanding of any numerical results. 
However, for the patterns A^, A 2 , B 3 and in Eqs. (9) and (10), we can reach the analytical results in 
a novel and simple way. The two zeros in these four patterns are located in the same row or colnmn of 
Mj^, namely, = 0 for a 7 ^ 6 . In this case, Eq. (4) can be recast into 


or equivalently 


Alt/al 

Ual 

+ 

TT 

Al^al 

As^aa 

Ubi 

+ 

^2^a2 TT 




/ ^l^al \ 

(Ki 

Ua2 


Ast^aa 
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^2Ua2 




-UaS, 

-Ub3-, 


(Ua3 

\Ut3 



(13) 


( 14 ) 
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Then, the orthogonality conditions of the unitary matrix U bring us another array of equations, which 
can be put in a form similar to Eq. (14), i.e., 



Ua2] 

U,2) 


fUk\ 

uh 

WJ 


(Ua3 

\Ut3 


(15) 


with c ^ a ^ b. If ~ ^a2^bi ^ which is generally satisfied for the allowed values of mixing 

parameters, the solutions to Eq. (14) and Eq. (15) should be identical. By identifying these two solutions, 
we obtain a set of new relations 


Al UclUa3 A, U*2Ua3 

As U*3Ual ’ Ag U*,U^2 ’ 


(16) 


from which one can extract two neutrino mass ratios C} and two Majorana CP-violating phases {p, a} 
in the same manner as in Eqs. (5) and ( 6 ). Obviously, the formulas in Eq. (16) are much simpler than 
those derived in previous works [301 EH [321 ng. 

As proved in Ref. [ 12 ], the predictions from and A 2 are related by a permutation symmetry 
corresponding to 6*23 90° — 6 * 23 , and likewise for Bg and B 4 . Thus, we apply Eq. (16) only to and 
B 4 , and discuss their phenomenological implications. 


• Pattern A 4 with = (Mj^)g^ = 0. Eor this pattern, we can specify the flavor indices in Eq. (16) 

as a = e, b = and c = t. Consequently, one arrives at 


A 3 

A 2 

A 3 


^*l^e3 I "13 

U^Uel C ?3 

^r2^e3 _ '^13 

U*r3Ue2 4, 


( ^ 12^23 ^\5 
V 02*^23 

( £l2£^gi<5 

V'*12^23 



e 


-2i(5 


e 


-2i(5 


(17) 


which are identical to those found in Ref. [12] up to different definitions of two Majorana CP-violating 
phases. To the leading order of S 43 , it has been found [30] 


^ ~ tan @23 tan ^42 sin ^13 , 

Q K. tan @23 cot 042 s™ ^13 ■ (18) 


With the help of Eqs. (11) and (18), one can get 0.98 < tan023 < 1-2, 0.64 < tan 042 < 0.68 and 
0.15 < sin 043 < 0.17 at the Icr level, implying that ^ < C < 1; namely, NO is favored. To the order 
of S 43 , both ^ and Q receive contributions from cos5, so we can compute 5 from Eq. ( 8 ), namely, 


cos 5 


tan 023 

tan 2042 sin 043 


/ sin 2042 tan 2042 p 
V 4 tan^ 023 sin^ 043 


It is now evident that 5 = 270° corresponds to 

^ ^ 4tan^ 023 sin ^043 
^ sin 2042 tan 2042 


(19) 


( 20 ) 


Given the best-fit values of = 7.50 x 10 ^ eV^ and = 2.58 x 10 ^ eV^ below Eq. (12), 

one obtains = 0.029. Using the best-fit values of three mixing angles, we immediately find that 


7 













the value on the right-hand side of Eq. (20) is 0.014, which is smaller than by a factor of two. 
Therefore, the best-fit values seem not to be consistent with the pattern A;^. 

As we have mentioned before, the results for A 2 can be derived from those for A^ by a replacement 
^23 90°-023- In this case, it is easy to verify that A 2 is compatible with NO as well. For 6 = 270°, 

tan^ 023 fhe right-hand side of Eq. (20) will be replaced by cot^ 023) implying a even smaller R^. 
Hence, we conclude that both A;^ and A 2 are consistent with NO, but they may run into problems if 
6 = 270° is experimentally conhrmed. However, if the current uncertainty on 6 is taken into acconnt, 
both and A 2 are not yet excluded. 


• Pattern B 4 with = 0. For this pattern, setting a = t, b = e and c = /r in Eq. (16), 

we arrive at 

^1 _ ^/j1^t3 _ C 23 '^12*^23 + 02'^23^13^ 

^^3^rl '®23 '®12'®23 “ '^12'^23'®13^'''''^ 

■^2 _ ^/j2^t3 _ ^23 02*^23 ~ '^12'^23^13^ (^21') 

'®23 02'523 + 

At the leading order, the neutrino mass ratios are given by ^ ~ ^ ~ cot^ 023, implying NO for 
023 > 45°. As for B 3 , the analytical results can be obtained by replacing with 90° — 023 in 
Eq. (21), which implies ^ k, tan^ 023 at the leading order. It is then straightforward to observe 
that only 10 is allowed for 023 > 45°. Therefore, we conclude that B 3 is now incompatible with the 
observed neutrino mass spectrnm. For B^, to the next-to-leading order, one can derive [T^] 


7^2 rrii 
m3 m3 


4cot^ 023 sin 0^3 
sin 20^2 sin 2023 


cos 5 , 


( 22 ) 


from which cos 5 < 0 is obtained because of m 2 > m^. Furthermore, the Dirac CP-violating phase 
is determined by 


cos 6 


sin 20^2 tan^ 023 
2 sin 0^3 tan 2023 


R„ 


(23) 


Using the best-fit values of three mixing angles and two neutrino mass-squared differences, one can 
get 5 ~ 269.6° or 90.4°, which is very close to the maximal CP-violating phase. The absolute neutrino 
mass scale is also fixed by neutrino mixing angles and mass-squared differences as 


m 3 = 


Amli 


1 — COt^ 023 


0.1 eV , m 2 ~ mi ^ m 3 cot^ 


^23 


0.087 eV 


(24) 


where 023 = 47° and Am^^ = 2.58 x 10“^ eV^ have been input. Obviously, neutrino masses are nearly 
degenerate, and their sum is S ~ 777 . 3(1 -|- 2cot^ 023 ) ~ 0.27 eV, which is already in tension with the 
current cosmological upper bound. In snmmary, the allowed region of the Dirac CP-violating phase 
is severely constrained, namely, only 6 ~ 90° or 270° is viable. If 023 ^ 45°, we have <5 —)■ 90° or 
270° and m^ —>■ 00 , so the upper bound on absolute neutrino masses forbids the maximal mixing 
angle 023 = 45° and maximal CP-violating phase 5 = 90° or 270° [35]. 


For other cases, Eq. (16) is no longer valid, and the analytical formulas for A 4 /A 3 and A 2 /A 3 become 
very complicated, as one can see from Refs. [301 ED Eg. To complete the list of two-zero textures that 
survive the latest neutrino oscillation data, we quote the analytical results for B 2 below 










• Pattern B 2 with = 0. For this pattern, we have 




'^12*^12*^2a(^'^23'^13 

' 5 i 2 '^ 12 C 23'®23 + (*^12 
■^12*^12*^23 (^'^ia'^ia 

'5i2*^12*^23'®23 + (*^12 


- C23C13) + ^23^13(^12^236'^ + ^12^236^'^) 

“ '5i2)'®23'^13®'^ + '5i2*^12*^23'®13(^ “I" ^23)^^'^ 

- C23C13) - ^23^13(^12^236*'^ + ^12^23^'^) 

“ '5i2)'®23'^13®''^ + '5i2*^12*^23'®13(^ + ^ 23)^^^*^ 


(25) 


from which one can obtain ^ Ri ^ « cot^ ^23 the leading order. Hence, B2 is compatible with NO 
for (923 > 45°, whereas ^ k, k, tan^ @23 is obtained for B^. The latter pattern is then excluded if 
the hints on NO and @23 > further confirmed by future neutrino oscillation data. At the 

next-to-leading order, the neutrino mass difference turns out to be 


7712 ^1 

7713 7773 


4sin0]^3 
sin 29^2 sin 2^33 


cos 6 , 


(26) 


indicating cos5 > 0 due to m 2 > mi- Moreover, the CP-violating phase is given by 


cos 6 « 


sin 2^12 

2 sin 0^3 tan 2^23 ‘ 


(27) 


Taking the best-fit values of neutrino mixing parameters, we get S ^ 270.3° or 89.7°. The absolute 
neutrino mass is the same as given in Eq. (24), and a nearly-degenerate neutrino mass spectrum is 
expected. To distinguish between B2 and B4, we have to measure the CP-violating phase as precisely 
as possible, e.g., an accuracy less than 1° is required if (5 ~ 270° is confirmed. 


Now we perform a thorough numerical analysis of all two-zero textures of by using the Icr ranges 
of neutrino mixing angles, Dirac CP-violating phase and neutrino mass-squared differences in Eqs. (11) 
and (12). The main strategy of our numerical analysis is as follows. (1) We randomly generate the values 
of { 6^2 ,013,023} and 6 in their la ranges, and then calculate two neutrino mass ratios {C,C} and i?^. (2) 
Several constraints are then chosen as criteria for whether a two-zero pattern is compatible with neutrino 
oscillation data. First of all, the fact of m 2 > m^ and NO, or equivalently ^ < C < 1; should be reproduced. 
Then, a pair of values for {AttiIi, AtttIi} are randomly generated and the ratio A77721 /A77731 is compared 
to Rj^ that has been computed from neutrino mixing angles and CP-violating phase in the first step. The 
difference between those two is required to be vanishingly small (e.g., < 10“^). The final numerical results 
are summarized in Fig. 1, and some comments are in order: 


1. Out of seven two-zero textures, only Ai 2 and B2 4 are found to be consistent with the neutrino 
oscillation data in Eqs. (11) and (12). In Fig. 1, the allowed regions of {012)^135^23} plotted with 
respect to those of 6 for A4 (red stars), A2 (blue squares), B2 (pink pluses), and B4 (green crosses), 
where the best-fit value 6 = 270° and the Icr interval 202° < 6 < 338° are represented by a dashed 
line and a shaded area, respectively. Different from Ref. [T5], the whole range of 6 is now considered. 

2. The patterns A4 and A2 are well separated in the parameter space of 6. In particular, the values 
of S for A4 are confined to a narrow area around 5 = 180°, i.e., the CP-conserving limit. Therefore, 
the sign of cos 5 (i.e., cos (5 < 0 for A4 and cos <5 > 0 for A2) is of crucial importance to discriminate 
between these two patterns. This is also the case for B2 and B4, as cos <5 > 0 for B2 and cos 6 < 0 
for B4. However, in the latter case, a very precise determination of 5 (e.g., less than 1°) is required. 

3. As one can observe from the lower plot of Fig. 1, a joint measurement of 023 and 5 is very useful to 
tell A2 apart from B24. For Ag, a nearly maximal CP-violating phase 5 ~ 90° or 270° is reachable 
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only if @23 takes a value close to the upper bound of its Icr range. Another way to distinguish 
between 2 and 63 4 is to discover neutrinoless double-beta decays, which never take place for the 
former cases because of a vanishing effective neutrino mass = | | = 0. But for B 2 and 

B 4 , we have = | \ ~ m 3 ~ 0.1 eV, which is likely to be probed in future experiments of 

neutrinoless double-beta decays [la Hi HU He]. 

The Dirac CP-violating phase 5 serves as an important discriminator for the two-zero textures of 
Although it seems difficult to further differentiate between B 2 and B 4 due to their similar implications at 
low energies, their high-energy phenomenologies in a realistic flavor model may be quite different. 

4 Symmetry Realization 

It is well known that the texture zeros in the Majorana neutrino mass matrix can always be realized 
by imposing Abelian flavor symmetries on the generic Lagrangian of neutrino mass models m- For 
instance, it has been demonstrated in Ref. [12] that cyclic symmetries can be used in the framework of 
type-II seesaw models [211 [22l [23l [Ml ESHShj to achieve all the viable two-zero textures of In general, 
however, it is a nontrivial task to derive two-zero textures by using non-Abelian discrete flavor symmetries. 
A predictive neutrino mass model with an flavor symmetry has been proposed in Ref. [H], where the 
patterns B^ and B 2 can be obtained. In this section, we recapitulate the essential idea to produce the 
two-zero texture of B 2 in Ref. [48], and then propose new scenarios to realize A^, A 2 and B 4 . 

In order to accommodate tiny neutrino masses, one can extend the SM by introducing three right- 
handed neutrinos (for i = 1,2,3) and a Majorana mass term for them. To further explain lepton 
flavor mixing, we can implement the non-Abelian discrete flavor symmetry |49l [50] and assign relevant 
fermion and scalar fields into suitable representations of the symmetry group. The assignments of fermion 
and scalar fields under the full symmetry SU(2 )l x U(1)y < 8 > A^ are 

4l ~ (2, -1) © 1 , 1 ', 1 " , ~ (1, - 2 ) © 3 , ~ ( 1 , 0 ) © 3 , ~ ( 2 , 1 ) © 3 , a ~ (3, - 2 ) © l" , (28) 

where the lepton doublets a = e, fi, r) are A^ singlets, three Higgs doublets form an A^ triplet, 

and one Higgs triplet A is an A^ singlet. Therefore, the gauge- and A 4 -invariant Lagrangian relevant for 
lepton masses and flavor mixing reads [38] 

+aj^ 

+ 2 ^^ (^r^r) ^ i^eh^^Th + + h-c- • (29) 

After the scalar fields acquire their vacuum expectation values (<hj) = v and (A) = u, the charged-lepton 
mass matrix and Dirac neutrino mass matrix Mj^ can be obtained from the first and second lines in 
Eq. (29), respectively. The effective neutrino mass matrix receives one contribution from the type-I 
seesaw formula —and the other from the triplet Higgs. In total, we have 



(^ 

0 



f al 

0 

0 ^ 

II 

0 

bA 

» 

-ii 

0 

0 

buCu 



0 

») 

Vo 

byC^ 

0 / 
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which can be identified as B 2 in Eq. (9) with four nonzero matrix elements a = —(Oj,n)^/M, b = a^u, 
c = b^u, and d = —b^c^v^/M. As shown in Ref. [38], if the Higgs triplet is assigned as 1', the pattern 
can be produced. In this model, however, it is impossible to realize any other viable two-zero textures. 

To illustrate how to produce the other two-zero textures, we consider the minimal type-(I-l-II) seesaw 
model m, where only one right-handed neutrino and one Higgs triplet A have been introduced to 
account for both tiny neutrino masses and baryon number asymmetry in our Universe. Now this minimal 
scenario is further extended, and the particle content and assignments under symmetry are 


Ul 






I'E. 


oR 


3 , A 


R 


1 , 




1 ' A 


(31) 


where two Higgs doublets if and (j) have been added as A 4 singlets. Except for a proper rearrangement 
of ^Th (^-®r)i'> charged-lepton sector remains unchanged and the mass 

matrix can be diagonalized just by rotating the right-handed charged-lepton fields in the flavor space, 
which has no impact on the lepton flavor mixing. The gauge- and A 4 -invariant Lagrangian relevant for 
neutrino masses becomes 


- ^ + O^^R^R + 9 ®A + ^-C. . (32) 


After the spontaneous symmetry breaking, the Dirac neutrino mass matrix turns out to be a 3 x 1 column 
vector and is given by M-^ = ( 0 , 019 , 6 ^)'^, where and bj^ = with = (ip) and = {(/)). 

The Majorana neutrino mass matrix is 



( ^ 

0 

«a\ ^ 

(0 

0 

0 ^ 

Mjj = u 

0 

bA 

» -F 

0 

al 



V^A 

0 

0/ 

VO 

Od^d 

bl / 


(33) 


which is exactly the pattern in Eq. (9) with a = a^u, b = b^u—a'^/M, c = —a-^b-^/M and d = —b'^/M. 
In a similar way, one can verify that the pattern A 2 can be achieved by assigning A as V under A 4 . Finally, 
to obtain B 4 , we need to make a modification of Eq. (31), i.e., 


4l ~ 1 ) ~ 1' ) 4l ~ 1'' ; ^aK ~3, $j~3, A~l, (34) 

leading to a Dirac neutrino mass matrix Mj-, = (ap, 6 ^, 0)"*" and thus the Majorana neutrino mass matrix 



(bA 

0 

0 ^ 

I 

“ M 

( «D 

Od&d 


M^ = u\ 

0 

0 

a A 


bl 



Vo 

a A 

0 ^ 

^ 0 

0 

0 / 


The phenomenological implications of these patterns for neutrino oscillation experiments have been ex¬ 
plored in previous sections. As we have seen, it is extremely difficult to distinguish between B 2 and B 4 
via neutrino oscillation experiments if the future oscillation data confirms a nearly-maximal CP-violating 
phase 5 ~ 270°. For instance, the accuracy of the determination of CP-violating phase is required to be be¬ 
low 1° for this purpose. In a specific model with flavor symmetries, however, new particles or interactions 
definitely cause additional observable effects that serve as a discriminator. 

We have not tried to write down the general scalar potential that is invariant under both gauge and 
flavor symmetries, and investigate whether the desired vacuum expectation values can be reached. This 
issne can be addressed in a similar way as in a number of works on applying the A 4 flavor symmetry to 
lepton flavor mixing [521 [53l I54j . Moreover, it is also interesting to examine if a successful leptogenesis 
works in our scenarios to explain the observed baryon number asymmetry in the Universe |55l 156] . These 
important questions deserve further dedicated studies, which are beyond the scope of our paper and left 
for fntnre work. 


II 







5 Summary 


In view of the latest results from atmospheric (Super-Kamiokande) and accelerator (T2K and NOz/A) 
neutrino oscillation experiments, which hint at the normal neutrino mass ordering < m 2 < m^, 
the maximal CP-violating phase 6 = 270° and the second octant ^23 > ^5°, we reexamine the two- 
zero textures of the Majorana neutrino mass matrix. By using the lo" ranges of neutrino oscillation 
parameters, we have found that only four patterns with two texture zeros (i.e., 2 ^ 2 , 4 ) pa-ss all 

current experimental constraints, and they can be further tested when the precision measurements of 
neutrino oscillation parameters are available. 

For the patterns 2 and Bg 4 , we have derived concise and useful formulas for neutrino mass ratios and 
Majorana CP-violating phases in Eq. (16), which have not yet been noticed in the literature. The simplicity 
of these formulas stems from the usage of unitarity conditions, and provides an excellent explanation for 
why the expressions in these four cases are much simpler than those in the cases of 2 and C. A detailed 
numerical analysis has also been carried out to check which two-zero texture is compatible with the latest 
experimental data. Our numerical results show that A^, A 2 , B 2 and B 4 survive the data at the la level, 
and precise determination of ^23 and 5 is of crucial importance to discriminate one pattern from another. 
Some distinct features should be mentioned: (1) The sign of cos S is positive for A 2 and B 2 , while negative 
for A 4 and B 4 ; (2) Only nearly-maximal CP-violating phase 5 ~ 90° or 270° is allowed for B 2 and B 4 , and 
the deviations are less than a few degrees; (3) The neutrino masses are nearly degenerate in the cases of B 2 
and B 4 , and substantial decay rates for neutrinoless double-beta decays are expected, while these decays 
do not occur in the cases of A^ and A 2 . Hence, in addition to oscillation experiments, the constraints 
on absolute neutrino masses and the observations of neutrinoless double-beta decays also help verify or 
disprove two-zero textures. As for model building, we have shown that the non-Abelian discrete flavor 
symmetry can be used to realize these four patterns in the type-(I-l-n) seesaw models. 

The ongoing and forthcoming neutrino oscillation experiments will unambiguously determine neutrino 
mass ordering and precisely measure the flavor mixing angles and the Dirac CP-violating phase. The tri¬ 
tium beta decays, neutrinoless double-beta decays and cosmological observations will probe the absolute 
neutrino mass scale and hopefully constrain the Majorana CP-violating phases. Any important experi¬ 
mental progress will give us new insights into the structure of lepton mass matrices, which can further 
help us explore the origin of neutrino masses and the true dynamics of flavor mixing and CP violation. 
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Figure 1: The allowed regions of neutrino mixing angles {^i 2 ) ^ 23 } versus the Dirac CP-violating phase 

S in the scenarios of (red stars), A 2 (blue squares), 63 (pink pluses), and B 4 (green crosses), where the 
shaded area corresponds to 202° < 6 < 338° and the vertical dotted line to 5 = 270°, and the lex ranges 
of neutrino mixing parameters in Eqs. (11) and (12) have been used. 
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